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Abstract 

We describe an elementary method to get non-asymptotic estimates 
for the moments of Hermitian random matrices whose elements are 
Gaussian independent random variables. We derive a system of recur- 
rent relations for the moments and the covariance terms and develop 
a triangular scheme to prove the recurrent estimates. The estimates 
we obtain are asymptotically exact in the sense that they give exact 
expressions for the first terms of l/7V-expansions of the moments and 
covariance terms. 

As the basic example, we consider the Gaussian Unitary Ensem- 
ble of random matrices (GUE). Immediate applications include the 
Gaussian Orthogonal Ensemble and the ensemble of Gaussian anti- 
symmetric Hermitian matrices. Finally we apply our method to the 
ensemble of N x N Gaussian Hermitian random matrices whose 
elements are zero outside of the band of width b. The other elements 
are taken from GUE; the matrix obtained is renormalized by 
We derive the estimates for the moments of and prove that the 

spectral norm H-ff^'^H remains bounded in the limit N, b — + oo when 
(log N) 3 / 2 /b^ 0. 



1 Introduction 



The moments of N x N Hermitian random matrices Hn are given by ex- 
pression 

Mf)= E {lTr(^}, 

where E{-} denotes the corresponding mathematical expectation. The asymp- 
totic behavior of Mj^ m the limit N —* oo is the source of numerous studies 
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and vast list of publications. One can observe three main directions of re- 
searches; we list and mark them with the references that are earliest in the 
field up to our knowledge. 

The first group of results is related with the limiting transition N — > oo 
when the numbers k are fixed. In this case the limiting values of , if 
they exist, determine the moments of the limiting spectral measure a of 
the ensemble {Hn}. This problem was considered first by E. Wigner [20]. 

Another asymptotic regime, when k goes to infinity at the same time as 
N does, is more informative and can be considered in two particular cases. In 
the first one k grows slowly and 1 <C k <C iV 7 for any 7 > 0. In particular, if 
k is of the order log iV or greater, the maximal eigenvalue of dominates 
in the asymptotic behavior of Then the exponential estimates of 

provide the asymptotic bounds for the probability of deviations of the 
spectral norm ||i2jv||- This observation due to U. Grenander has originated 
a series of deep results started by S. Geman [1, 7, 9]. 

The second asymptotic regime is related to the limit when k = 0(iV 7 ) 
with 7 > 0. The main subject here is to determine the critical exponent 7 
such that the same estimates for M%jp as in the previous case remain valid 
for all 7 < 7 and fail otherwise [18]. This allows one to conclude about the 
order of the mean distance between eigenvalues at the border of the support 
of the limiting spectral density da [4, 19]. 

In present article we describe a method to get the estimates of M^p that 
are valid for all values of iV and k such that k < CN^ with some constant 
C. The estimates of this type are called non-asymptotic. However, they 
remain valid in the limit ./V — ► 00 and in this case they belong to the second 
asymptotic regime. 

As the basic example, we consider the Gaussian Unitary (Invariant) En- 
semble of random matrices that is usually abbreviated as GUE. In section 
2 we describe our method and prove the main results for GUE. Immedi- 
ate applications of our method include the Gaussian Orthogonal (Invariant) 
Ensemble of random matrices (GOE) and the Gaussian anti-symmetric (or 
skew-symmetric) Hermitian random matrices with independent elements. 
Detailed description of these ensembles is given in monograph [16]. In sec- 
tion 3 we present the non-asymptotic estimates for the corresponding mo- 
ments. 

Our approach is elementary. We use the integration by parts formula 
and generating functions technique only. We do not employ such a powerful 
methods like the orthogonal polynomials technique commonly applied to 
unitary and orthogonally invariant random matrix ensembles. This allows 
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us to consider more general ensembles of random matrices than GUE and 
GOE. One of the possible developments is given by the study of the ensemble 
of Hermitian band random matrices H^ N ' b \ The matrix elements of H^ N,b ^ 
within the band of the width b along the principal diagonal coincide with 
those of GUE. Outside of this band they are equal to zero; the matrix 
obtained is normalized by fe -1 / 2 . In section 4 we prove non- asymptotic 
estimates for the moments of H^ N ' b ^ . These estimates allow us to conclude 
about the asymptotic behavior of the spectral norm \\H( N ^\\ in the limit 
b, N -»■ oo. 

In section 5 we collect auxiliary computations and formulas used. 



1.1 GUE, recurrent relations and semi-circle law 

GUE is determined by the probability distribution over the set of Hermitian 
matrices {H^} with the density proportional to 

exp{-27V TiH^}. (1.1) 

The odd moments of Hn are zero and the even ones M^jp verify the following 
remarkable recurrent relation discovered by Harer and Don Zagier [11] 



(N) _ 2k -1 {N) 2k -1 2fc-3 fc(fc-l) (at) 
M2fc " 2k + 2 M2k ~ 2 + 2k + 2 2k 4AT2 M ™~*> 



where Mq N ^ = 1 and = 1/4. It follows from (1.2) that the moments 

M^\k = 0,1,... converge as N — > oo to the limiting mu determined by 
relations 

2k — 1 

m k = 2k + 2 m k-i, m = 1. (1.3) 

The limiting moments {m^, k > 0} are proportional to the Catalan numbers 
Ck- 

1 1 (2k\ 1 ^ . . 

mk = ¥WTT){k)=¥ Ck (L4) 

and therefore verify the following recurrent relation 

j fc-i 

m k = j m k-i-j rrij, A; = 1,2, ... (1.5) 

j=0 

with obvious initial condition mo = 1. 
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In random matrix theory, equality (1.5) was observed for the first time 
by E. Wigner [20]. Relation (1.5) implies that the generating function of the 
moments 

oo 

/(r) = J2 m k-r k 
verifies quadratic equation r/ 2 (r) — 4/(r) + 4 = and is given by 



Using (1.6), Wigner has shown that the measure a w determined by the 
moments = J \ 2k da w (X) has the density of the semicircle form 



if|A|<l, (17) 
^ (A) -vr\0, if|A|>l. (L7) 

The statement that the moments converge to for I = 2k and to 

for I = 2k + 1 is known as the Wigner semicircle law. 

In present paper we show that the generating function /(r) together 
with its derivatives represents a very convenient tool when estimating the 
moments M%p . Everywhere below, we use denotation [•]& for the fc-th 
coefficient of the corresponding development, so [f{r)\k = rn*k- 

1.2 Estimates for the moments of GUE 

Using relations (1.2) and (1.3), one can easily prove by induction the esti- 
mate 

^JP<(l + ^y rn k . (1.8) 

Indeed, let us assume inequalities M 2 ( f } < (1 + I 2 / (gN 2 )) 2l mi with some 
g > to hold for all values of I such that 1 < I < k — 1. Let us show that 
this is also true for / = k provided g < 8. 

Regarding the right-hand side of (1.2) and replacing M^l 2 an d ^2k-4 
by corresponding estimates with I = k — 1 and I = k — 2, respectively, we 
bound the right-hand side of (1.2) by the sum of 

2k- 1 
2k + 2 
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and 



Here we have used identity (1.3). Comparing the expression obtained with 
the right-hand side of (1.8), we see that the following inequality 





is sufficient for (1.8) to be true. Expanding the powers, we see that the 
condition g < 8 is sufficient to have (1.8) valid for all values of k and N. 

Estimates (1.8) are valid for all values of k and ./V without any restric- 
tion. They allow one to estimate the probability of deviations of the largest 
eigenvalue of Hn (see, for example [14, 15] and references therein). Then 
one can study the asymptotic behavior of the maximal eigenvalues and also 
conclude about spectral scales at the borders of the support of a' w (see [18]). 

It should be noted that relations (1.2) are obtained in [11] with the help 
of the orthogonal polynomials technique (see [10] and [15] for the simplified 
derivation). There are several more random matrix ensembles (see [15] for 
the references) whose moments verify recurrent relations of the type (1.2). 
But relations of the type (1.2) are rather exceptional than typical. Even in 
the case of GOE, it is not known whether relations of the type (1.2) exist. 
As a result, no simple derivation of the estimates of the form (1.8) for GOE 
has been reported. 

We develop one more approach to prove non-asymptotic estimates of the 
type (1.8). Instead of relations (1.2), we use the system of recurrent relations 
(1.5) that is of more general character than (1.2). Regarding various random 
matrix ensembles, one can observe that the limiting moments verify either 
(1.5) by itself or one or another system recurrent relations generalizing (1.5) 
(see for instance, section 5 of [3], where the first elements of the present 
approach were presented). 

We derive a system of recurrent relations for the moments M^P that 
have (1.5) as the limiting form. These relations for M^jp involve correspond- 
ing covariance terms. Using the generation functions technique, we find the 
form of estimates and use the triangle scheme of the recurrent estimates to 
prove the bounds for moments and covariance terms. The final result can 
be written as 




(1.9) 
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with some a > 1/12. The estimates obtained are valid in the domain k 3 < 
xN 2 with some constant x, le. not for all values of k and N, as (1.8) does. 
But in this region our estimates are more precise than those of (1.8). If 
k 3 <C N 2 , our estimates provide exact expressions for 1/iV-corrections for 
the moments M^P . 

1.3 Band random matrices and the semi-circle law 

Hermitian band random matrices H^ N,b ^ can be obtained from GUE matrices 
by erasing all elements outside of the band of width 6 along the principal 
diagonal and by renormalizing the matrix obtained by the factor 6 -1 / 2 . It 
appears that the limiting values of the moments 

M^= E {lTr(^)) 2fe } 

crucially depend of the ratio between 6 and ./V when N — > oo (see [5, 13, 17]). 

If b/N — ► 1 as N — > oo, then M^'^ — > rrik and the semicircle law is 
valid in this case. If b/N — ► c and < c < 1, then the limiting values of 
M^'^ differ from m^. Finally, if 1 <C b <C AT", then the semicircle law is 
valid again. 

The last asymptotic regime of (relatively) narrow band width attracts a 
special interest of researchers. In this case the spectral properties of band 
random matrices exhibit a transition from one type to another. The first 
one is characterized by GUE matrices and the second is given by spectral 
properties of Jacobi random matrices, i.e. the discrete analog of the random 
Schrodinger operator with 6 = 3 (see [6, 8] for the results and references). 
It is shown that the value b' = \fN is critical with respect to this transition 
[6, 8, 12]. 

In present paper we derive the estimates for M^'^ that have the same 
form as the estimates for GUE with N replaced by b. This can be viewed as 
an evidence to the fact that the asymptotic behavior of the eigenvalues of 
jj(N,b) a £ border of the semi-circle density is similar to that of matrices 
of the size 6x6. The estimates we obtain show that the value 6' = \f~N~ 
does not play any particular role with respect to the asymptotic behavior 
of the spectral norm WH^^W. We show that if 6 > (logiV) 3 / 2 , then the 
spectral norm converges with probability 1 when iV — > oo to the edge of 
the corresponding semicircle density. Up to our knowledge, this is the first 
result on the upper bound of the spectral norm of band random matrices. 
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2 Gaussian Hermitian Ensembles 



(H™) xy = -=h xy , x,y = l,...,N (2.3) 



Let us consider the family of complex random variables 

, _ jV xy + iW xy , iix<y, 
xy ~ \V yx - iW yx , ifx>y, 

where {V xy ,W xy , 1 < x < y < N} are real jointly independent random 
variables that have normal (Gaussian) distribution with the properties 

BV xy = BW xy = 0, (2.2a) 

and 

EV* y = (l + 5 xy )±±lL, EW* y = (l-6 xy )^>, (2.26) 

where 5 xy is the Kronecker (5-symbol and r/ £ [—1, 1]. Then we obtain the 
family of Gaussian ensembles of N x N Hermitian random matrices of the 
form 

frh - -L 

that generalizes the Gaussian Unitary Ensemble (1.1). Indeed, it is easy to 
see that {H$} coincides with the GUE, while and {H^} repro- 

duce the GOE and Hermitian skew-symmetric Gaussian matrices. In [16], 
the last ensemble is referred to as the Hermitian anti-symmetric one; below 
we follow this terminology. The present section is devoted to the results 
for GUE and their proofs. Two other ensembles will be considered in the 
section 3. 

2.1 Main results for GUE and the scheme of the proof 

Let us consider the moments of GUE matrices. We prove a little more 

precise estimate than (1.9). 

Theorem 2.1 

Given any constant a > 1/12, there exists x > such that the estimate 

M^<(l + a k -^-y k (2.4) 

holds for all values of k,N under condition that k 3 /N 2 < \- 

Remark. Using relation (1.2), one can prove (2.4) under condition that 

a >^- x (2.5) 
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This relation shows that Theorem 2.1 gives the correct lower bound for a. 
In our proof we get relations between x an d a more complicated than (2.5), 
but they are of the same character as (2.5). It follows from (2.5) that the 
closer a to 1/12 is, the smaller % has to be chosen and vice versa. Indeed, 
the following proposition shows that the estimate (2.4) is asymptotically- 
exact. 

Theorem 2.2 

Given k fixed, the following asymptotic expansion holds, 

M 2k ] =™fc + ^™fc 2) +0(N- A ), asN^oo, (2.6a) 



where 

i(2 ) _ fc(fc-l)(fc + l) 

12 

// k — ► oo and % = k 3 /N 2 — > 0, i/ien relation (2.6a) remains true with 
0(N~ 4 ) replaced by o(x). 



Remark. It follows from (1.2) that the sequence {m^ 2 -*, fc > 1} is determined 
by recurrent relation 

ro< 2 > - ^± ■ m< 2 > fc(fc " 1) mh fc-12 
with obvious initial condition = 0. It is easy to check that (2.6b) is in 

(2) 

complete agreement with this recurrent relation for m k . 

Let us explain the role of recurrent relations (1.5) in the proof of Theorem 
2.1. To do this, let us consider the normalized trace L a = Tr H a 

1 N 

x,s=l 

and compute the last mathematical expectation. Here and below we omit 
subscripts and superscripts iV when no confusion can arise. Applying the 
integration by parts formula (see section 5 for details), we obtain equality 

E{M = -^E{W 3 L 3 }. (2.7) 

3=0 
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Introducing the centered random variables L° = Lj — ELj, we can write 
that 

E {L ai L a2 } = E{L ai }E{L a2 } + E{L° ai L° a2 }. 
Taking into account that EL 2 k+i = 0, we deduce from (2.7) relation 

<> = \ £ M%1 2 _ 2) M«f > 4- \ fl£3, (2.8) 

3=0 

where we denoted 

ai+a,2=2k—2 

Obviously, the last summation runs over a« > 0. Comparing (2.8) with (1.5), 
we see that the problem is to estimate the covariance terms D^ 2 \ Here and 
below we omit superscripts N when no confusion can arise. 

In what follows, we prove that under conditions of Theorem 2.1, 

i«r*i < (2.9) 

with some constant c. Inequality (2.9) represents the main technical result 
of this paper. It is proved in the next subsection. With (2.9) in hands, we 
can use relation (2.8) to show that (2.4) holds. 

Now let us explain the use of the generating function /(r) (1.6). Regard- 
ing the right-hand side of (2.4), one can observe that the third derivative 
of /(t) could be useful in computations because of the equality [/ w (t)] = 
(k + 3)(k + 2)(k + l)mfc + 3. Indeed, more accurate computations (see identity 
(5.12) of section 5) show that the function 

A t 2 So 
f(r) + W2 jj—^ = * N { T ) with A = - (2.10a) 

is a very good candidate to generate the estimating expressions. This is 
not by a simple coincidence or an artificial choice. Later we will see that 
the form of <&n{t) is in certain sense optimal. It is dictated by the iteration 
scheme we use to get l/7V-corrections for the moments and covariance terms 
(see subsection 2.5, the proof of Theorem 2.2). 
Let us now show how (2.9) implies the estimate 

M^<[^ N (r)] k . (2.106) 
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Assuming that this estimates and (2.9) are valid for all the terms of the 
right-hand side of (2.8), we can estimate it with the help of inequalities 



fc-i 



£ M 2fc _ 2 _ 2j M 2j + - A \Dfl 2 \ < - *%(r) 



+ 



3=0 



k-i 4N 2 



(1-r) 



fc-2 



Denoting by 6(fc; iV) the terms of the order 0(N 4 ), we can write that 

A 



rf 2 (r) + r 3 /(r) 



2 JV 2 (1-t) 5 / 2 



+ 0(A;;JV). 



Rewriting (1.6) and quadratic equation for /(t) in convenient forms 

i-i/t^f, (2.ii) 



^»=/M-l and T >« 



4 ' v ' 2 
we transform the expression in the brackets: 



f(r) + 



A 



N 2 (1 - r) 5 /2 A^ 2 (1 - r) 2 



A 



(1-r) 



Remembering that [<E>7v( r )]fc reproduces the expression to estimate 
we conclude that (2.10) is valid provided 



A 

iV 2 



{l-rf 



> 



4AT2 



{l-rf 



cjk - 1) 
4iV 2 



(2.12) 



This requires inequality A > c/4. 

The final comment is related to the role of the terms Q(k; N). They are 
of the form 



&{k;N) 



A 2 

4AT4 



(1-r) 5 



< 



A 2 k A 

AT4 • 



If one wants these terms not to violate inequality (2.12) involving terms of 
the form k/N 2 , one has to set the ratio k s /N 2 = x sufficiently small. This 
explains the last condition of Theorem 2.1. 

It should be noted that the same comments concern the proof of the 
estimate of covariance terms (2.9), where the recurrent relations, generating 
functions and terms of the type x appear. In the proofs, we constantly use 
relations (2.11). 



10 



2.2 Main technical result 

In this subsection we prove the estimates of the covariance terms of the type 
D^k = ^^{^1^2}' ^he mam idea is that these terms are determined by 
a system of recurrent relations similar to (2.8). These relations involve the 
terms of more complicated structure than . The variables we study are 
defined as 



D 



(?) 

2k 



E D it,a q = E e{z4l° 2 ...l°J, q>2. 



aj+...+a„=2k 



ai+...+a„=2k 



Here and everywhere below, we assume that the summation runs over all 
positive integers a« > 0. 

Our main technical result is given by the following statement. 



Proposition 2.1. 

Given A > 1/16, there exists x > such that estimate (2.10) holds for all 
values of 1 < k < kQ, where ko verifies condition 



h 3 
N 2 



<X- 



Also there exists C 



such that inequalities 



1 2A 

— <C<max{ — ,4!} 



,(3*)! 



and 



\ Di 2k S) \<C- N2s 



i n (2»+l), < r ( 3s + 3 ) ! 
\^2k I — N 2s+2 



(1 - r) 2s 



(1 _ T )2s+5/2 



are irue /or all k,s such that 

2k + q< 2k 
with q = 2s and q = 2s + 1, respectively. 



(2.13) 
(2.14) 

(2.15a) 
(2.156) 

(2.16) 



Remark. The form of estimates (2.15) is dictated by the structure of the 
recurrent relations we derive below. The bounds for the constants A and C 
and of the form factorial terms of (2.15) are explained in subsection 2.4. 
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We prove Proposition 2.1 in the next subsection on the base of recurrent 
relations for that we derive now. Let us use identity for centered random 
variables E,{X°Y°} = E{XY°} and consider equality 



E{K 1 LZ 2 ...LZ q }=E{L ai [L^---L° aq }°} = ^ £ E {h^H^ 1 ^ ■ ■ ■ K, 

x,s=l 

(2.17) 

We apply to the last expression the integration by parts formula (5.1) and 
obtain equality 

Di at.,a q = li:v{La 1 -2- j L j [L° a2 ...L° aq }°} 



3=0 

1 r 

-f 



1 9 

M2 £ E { L a2 • • • L °ai-i °* ^i+a^K.^ . . . L° q } , (2.18) 



AN 2 

i=2 

with the help of formulas (5.7) and (5.8), respectively. The detailed deriva- 
tion of (2.18) is presented in subsection 5.2. 

Let us consider the first term from the right-hand side of (2.18). We can 
rewrite it in terms of variables D with the help of the following identity 

E{L!L 2 Q°} = E{Li}E{L§Q}+E{L 2 }E{L?Q}+E{L?L§Q}-E{L?L§}E{Q}, 

where Q = L° 2 ■ ■ ■ L° q . Regarding the last term of (2.18), we use (2.17) and 
obtain relation 

-. a\— 2 -i oi— 2 

n(q) _ _ m n(«) i _ m n( ? ) 

±J a 1 ,...,a q - a 2^ 1V1 3 ly a 1 -2-j,a 2 ,...,a q f , L -^-3^3,02 ,...,a q 

3=0 3=0 

+i V - - V D {2) r)(i^ 1 ) 

1 4 Z-^ j,a\-2-j,a2,—,a q a Z^i j,ai-2-j a 2 ,...,a q 
3=0 3=0 

1 9 1 9 ( i) 

^4JV2 £ a i^ai+a i -2-D a 2,...,a l _i,a I+1 ,...,a 9 + 4 jY2 £ a *^ ) a 2 ,...,a i _ 1 ,Oi+oi -2,o i+ i ,...,a 9 • 
i=2 i=2 

(2.19) 

Taking into account that M^_ 1 = 0, it is easy to deduce from (2.19) by 
induction on k that 

D ai,...,a q = whenever ai + . . . + a q = 2k + 1. 



12 



Let us introduce variables 

D$ = E Dil] 



a q 

cii + ...+a q =2k 



Using the positivity of M 2 j, we derive from (2.19) the second main relation 



2 ~2fc-2-2^"<y ' ^ N 2 Z 

where 1 < k, 2 < g < 2k. When regarding two last terms of (2.19), we have 
used obvious equality 



-F 



o'-2- 



«i+a2=a' \o2=l / 

Using this relation with F replaced by M and a' = 2j + 2, we obtain that 

V „ m I n(9) l - V (2j + 2)(2j + i) 

2-, ^■ M ai+o i -2|-L'iV ) ... > a i _ 1)ai+ll ... ) o,l ~ 2_, D 2k-2-2j g M ^3 ■ 

ai + ...+a q =2k j=0 

Also we can write that 



i+ai—2,a i+1 ,...,a q I 

a\ + ...+a q =2k 



V a\D {q - 1] 



E i<lix E (^-«i+2)< E i^fcVJx^-^ 3 

6 2 +...+b q =2fc-2 l<ai<6i+l ai + ...+a 9 _i=2fc-2 

and get the last term of (2.20). 

The upper bounds of sums in (2.20) are written under agreement that 

Dyj = whenever q > 2k. Also we note that the form of inequalities (2.20) 
is slightly different when we consider particular values of q and k. Indeed, 
some terms are missing when the left-hand side is D^u, D^? , D^t ^> 

— (2k— 2) — (a) 

and D 2k . However, the agreement that Dj^ = whenever q > 2k and 

that = and = 5k,o make (2.20) valid in these cases. 
Obviously, we have that 



1 1 - 

M 2k < - J2 M 2jfc _ 2 _ 2j M 2j + -Dfl 2 . (2.21) 

3=0 
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2.3 Recurrent relations and estimates 

To estimate M and D^ q \ we introduce auxiliary numbers > and 

R^k^ — determined by a system of two recurrent relations induced by 
(2.20) and (2.21). This system is given by the following equalities (we omit 
superscripts N) 

B k = \{B *B) k _ 1 + \r^v (2-22) 

and 

+ I(ij(.-)»i J W) fc _ i + ^4r i 1) , (2.23) 
considered in the domain 

A = {(k,q) : fc> 1, 2 < g< 2fc} 

with denotation 

(2fc + 2)(2fc + l) 
^fc 2 fc 

and the convolutions as follows 

fc-i 

(^*^) fe -i = E^-i-^- 

3=0 

The initial values for (2.22)-(2.23) coincide with those of M and D: 

bP = 1, R?< N) - — 



AN 2 ' 



Let us note that one can consider relations (2.22) and (2.23) for all integers 
k and q with obvious agreement that outside of A the values of R are zero 
except the origin Rp N "> = l. The system (2.22)-(2.23) plays a fundamental 
role in our method the proof of Proposition 2.1. This proof is composed of 
the following three statements. 

Lemma 2.1. 

Given fixed N, the family of numbers {B^,R^\ (k,q) £ A} exist; it is 
uniquely determined by the system of relations (2. 22) -(2. 23). 

Lemma 2.2. 
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Inequalities 



< B™ 



and D^ N) <Rf N) 



(2.24) 



/ioW /or a// N and (k, q) 6 A. 



Lemma 2.3. 

Under conditions of Proposition 2.1, the numbers and are estimated 
by the right-hand sides of inequalities (2.10), (2.14) and (2.15), respectively; 
that is 



B[ N) < 



f{T)+AN- 2 T 2 {l-T) 



-5/2 



and 



C{3s)!N- 2s [r(l -t)~ 2s \ 
C(3s + 3)\N- 



-2s-2 



k ' 



ifq = 2s; 



r(l - r )-(4s+5)/2 > if q = 2s + l. 

. k 



(2.25) 



(2.26) 



Lemma 2.3 represents the main technical result concerning the system 
(2.22)-(2.23). Lemmas 2.1 and 2.2 looks like a simple consequence of the 
recurrent procedure applied to relations (2.22)-(2.23) and (2.20)-(2.21), re- 
spectively. However, the form of recurrent relations (2.22)-(2.23) is not usual 
because relations for B involve the values of R and vice-versa. The ordinary 
scheme of recurrence has to be modified. This modification is described in 
the next subsection. Lemma 2.3 is also proved on the base of this modified 
scheme of recurrence. 



2.3.1 The triangular scheme of recurrent estimates 

Let us show on the example of Lemma 3 that the ordinary scheme of re- 
current estimates can be applied to the system (2.20)-(2.23). Under the 
ordinary scheme we mean the following reasoning. Assume that the esti- 
mates we need are valid for the terms entering the right-hand side of the 
inequalities derived. Apply these estimates to all terms there and show that 
the sum of the expressions obtained is smaller than that we assume for the 
terms of the left-hand side; check the estimates of the initial terms. Then all 
estimates we need are true. Let consider the plane of integers (k, q) assume 
that estimates (2.26) are valid for all variables R with (k, q) lying inside of 
the triangle domain A(m),m > 3 

A(m) = {(k, q): 1 < k, 2 < q < 2k, k + q < m} 
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and that estimates (2.25) are valid for all variables B\ with 1 < / < m — 2. 

Then we proceed to complete the next line k + q = m + 1 step by step 
starting from the top point T{m + 1) of A(m + 1) and ending at the bottom 
point (m — 1, 2) of this side line. This means that on each step, we assume 
estimates (2.25) and (2.26) valid for all terms entering the right-hand sides 
of relations (2.23) and show that the same estimate is valid for the term 
standing at the left hand side of (2.23). 

Once the bottom point (m— 1, 2) achieved, we turn to relation (2.22) and 
prove that estimate (2.25) is valid for B m . Again, this is done by assuming 
that all terms entering the right-hand side of (2.22) verify estimates (2.25) 
and (2.26) with q = 2, and showing that the expression obtained is bounded 
by the right-hand side of (2.25). This completes the triangular scheme of 
recurrent estimates. 

It is easy to see that the reasoning described above proves, with obvious 
changes, Lemmas 2.1 and 2.2. 



2.3.2 Estimates for B 



Assuming that the terms standing in the right-hand side of (2.22) are esti- 
mates (2.25) and (2.26) with s = 1, we can write inequality 



\(B*B) k . 1 + \R^. 1 



< 



rf 2 (r) | A 



r 3 /(r) 



A 2 



3C 



N 2 2(1 - r) 5 / 2 + 4iV 4 (1 - r) 5 fc + 2N 2 (1 - r) 2 fc " 

Taking into account relations (2.11), we transform the first bracket of 
(2.27) into expression 



f(r) + 



A 



N 2 (1 - r) 5 /2 



J k 



A_ 

N 2 



(1 



+ 



J k 



A 2 
4iV 4 



(1 



J k 



Here, the first term reproduces expression [$jv(T)]fc ; the second term is 
negative and this allows us to show that the estimate wanted is true. Then 
we see that estimate B^ < [$Jv(T)]fc 1S true whenever inequality 



A 



(l-T)S 



> 



3C 



{l-rf 



+ 



A 2 
AN 2 



(l-T)E 



(2.28) 



holds. This is equivalent to the condition 



A > ^ + A " (fc-4)(fc-3)(fc-2) 
" 2 ^ AN 2 4! 
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Remembering that k 3 < xN 2 , we see that the estimate (2.25) of is true 
provided 

A > — + ^2L. (2.29) 

~ 2 96 ' 



2.3.3 Estimates for R^ 

Let us rewrite (2.8) with q = 2s, s > 2, k > 1 in the form 

nt > = 1 (^-) * B ) k _ i + *^ * i + x + y + z, (2.30) 

where we denoted 

X = iie^\ F = l(i ? (—)*^)) fe _ i , Z = ^£r 1 ). (2.31) 

The first term in the right-hand side of (2.30) admits the following esti- 
mate 



C(3s)\ 

N 2s 



+ 



A 



r 2 f(r) 

2(1 - t) 2s ' 2iV 2 (1 - r) 2s + 5 /2 



Using (2.11), we transform the last expression to the form 



C(3a)! 



+ 



4 



(l-r) 2s (l-r) 2s -V2 2A^ 2 (1 - r ) 2s + 5 / 2 



(2.32) 



The first term reproduces the expression we need to estimate R^ s ^ ■ 

Let us consider the second terms of the right-hand side of (2.30). As- 
suming (2.25) and using identities of subsection 5.1, it is not hard to show 
that 

1 18A 1 1 . . 

+ — — . (2.33) 



_(l_ T )3/2 ^2 (l_ r )9/2 

Indeed, it follows from (5.11) that 

(2fc + 2)(2fc + l) r ,, ,. {2k + 2){2k + l) 
2 U( r )Jfc = 7, m k = 

Next, identity (5.12) implies relation 



(1 - r) 3 /2 



(2fc + 2)(2fc + l) 



(1 - r) 5 /2 



(2fc + 2)(2fc + l) 2k(2k- l)(2fc + 1) 
2 3! 
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Now, regarding (5.9) with r = 4, it is easy to see that 

(2k + l)(2k + 2)(2k + 3)(2ife + 5)(2fc + 7) 



(1 - r) 9 / 2 



5-6-7 



m fe . 



Then (2.33) follows. 

Returning to the right-hand side of (2.30), we can write with the help of 
(2.33) inequality 



(R^ S - 2 U B") k i < 



C(3s - 3)! 
4jy2s-2 



+ 



ISA 



(1 _ T )2s-l/2 AT2 (1 _ T )2s+5/2 



(2*34) 

Here and below we use relation [t-^t)]*, < [g(T)]k valid for the generating 
functions under consideration. Let us stress that (2.34) remains valid in the 
case of s = 1 with C replaced by 1. 

Let us turn to (2.31). We estimate the sum of X and Y by 



X + Y < 



C(l + C)(3s + 3)! 



For the last term of (2.31) we can write inequality 

Ck 2 (3s + 1)! 



(1 _ r )2s+5/2 



Z < 



(1 _ r )2s+l/2 



(2.35) 



(2.36) 



Comparing the second term of (2.32) with the sum of the last term of 
(2.32) and the right-hand sides of (2.34), (2.35), and (2.37), we arrive at the 
following inequality to hold 



c > (2s-l)(3s-3)! _ 8s± + C(l - M k 2 (3s + 1) 



r(l-r) 



-2s-l/2 



(3s)\ 4 3iV 2 



[t(1- 



-2s+l 



(1 + C)(3a + 3)! + 18A(3s - 2)! + 2A(3s) 



r (l _ r )-2-5/2- 



4iV 2 (3s)! 
Using identity (5.10), we see that 

T (l _ r )-2.-l/2 



[ T (l _ r )-2 S +l/2] 



2.37) 



Similarly 



[t(1-t)-W\ 
t(1 - t)- 2s - 5 / 2 



2k + As - 2 < 4A; 



4s - 1 



4s - 1 



< 



(4fc ) 3 



^(I-t)- 2 ^ 1 /^^ - (4s-l)(4s + l)(4s + 3)' 
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Inserting these inequalities into (2.37), maximizing the expressions obtained 
with respect to s, and using (2.13), we get the following sufficient condition 



C > ^ + °^ — ^ + 2 X (1 + 10C(1 + C) + 2AC) . 



(2.38) 



2.3.4 Estimates for J R( 2s+1 ) 

Let us turn to the case q = 2s + 1 and rewrite (2.8) in the form 

(2.39) 

where 

X, = \RtT\ Yl = \ (BP-* * RM)^ , * = (2.40) 
Regarding the first term of (2.39), we can write inequality 

r 2 /(r) 



2\ Jk- N 2s + 2 



+ 



At 



C(3s + 3)! 

]\[2s+2 



2(l- r )2 S +5/2 2iV 2 (l - T ) 2s + 5 

At 



+ 



(l_ r )2s+5/2 ( 1 _ T )2s+2 1 2N 2 (1 - r) 2s + 5 



. (2.41) 



The first term of the right-hand side of (2.41) reproduces the expression 
needed to estimate R^ s+1 \ 

Let us consider the second term of (2.39). It is estimated as follows: 



(R^-V* B") ^ < 



Cs(3s)\ 



_ t-\2s+2 



(1-T 



9ACs(3s)} 

N2s+4 



2iV 2 V Jk-i ~ 2N 2s + 2 

Regarding two first terms of (2.40), we can write that 
C(3s + 3)! + 6C 2 (3s)! 



(1-r) 



2s+5 



X 1 +Y 1 < 



and 



Zi < 



4]\[2s+2 

Ck 2 s(3s)\ 

jy2s+2 



(1 _ r )2s+2 



(1-r) 



2.s 



Comparing the negative term of (2.41) with the sum of the last term of 
(2.41) and the estimates for the terms of (2.40), we obtain inequality 



c(% 



s(3s)\ 



> 



3C 2 (3s)! k 2 s(3s)\ [r(l-r) 



-2sl 



+ 



4 2(3s + 3)\J ~ 2{3s + 3)! (3s + 3)! [r(l - r)^- 2 ] fe 
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AC / 18,(3,)! N [r{l-r)^-% 

+ 2N 2 { l+ {3s + 3)] ) [r{1 _ T) -2s-2 ]k ^ 

Equality (5.13) implies that 

[t(1-t)- 2 % 2s(2s + 1) 



and that 



[ T (l _ r )-2 S -2] fc l + 2s)(A; + 2s) 



[r(l - r)" 2 - 8 ]* 8fc 3 



[r(l - T)' 2s ' 2 } k ~ (2s + 2)(2s + 3)(2s + 4) ' 

Inserting these two relations into (2.42) and maximizing expressions with 
respect to s, we obtain, after elementary transformations, the following suf- 
ficient condition 

4' 

C < i T7-- ( 2 - 43 ) 

" 1 + AA X 



2.4 Proof of Theorem 2.1 

Let us repeat that inequalities (2.29), (2.38), and (2.43) represent sufficient 
conditions for recurrent estimates (2.25) and (2.26) to be true. Let A > 1/16. 
Then for any constant C < 4! verifying condition 

1 r , 2A 
24 <C< 1P 

there exists such x > that (2.38) is true. Indeed, it is sufficient to take 
X < x'i where %' is such that 

1 23 

2 X 'K < min{C ,—C\, 

i 24 24 ; 

with K = 1 + 10(7(1 + C) + 2AC. Also there exists x" such that (cf. (2.29)) 

A > 3C/2 + A 2 X ". 

The choice of X — mm {x') x"} makes (2.29) and (2.38) true. Condition 
(2.43) is obviously verified. Thus, conditions (2.13), (2.14), and A > 1/16 
of Proposition 2.1 are sufficient for (2.29), (2.38), and (2.43) to hold. This 
completes the proof of Lemma 2.3. 

Lemma 2.2 together with Lemma 2.3 implies estimates (2.10) and (2.15). 
Then Proposition 2.1 follows. The statement of Theorem 1.1 is a simple 
consequence of the estimate (2.10) and Proposition 2.1. 
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We complete this subsection with the discussion of the form of estimates 
(2.26) and constants A and C. First let us note that the upper bound 4! 
for C imposed by (2.14) represents a technical restriction; it can be avoided, 
for example, by modifying estimates (2.26) for i?( 2s ) and R( 2s+1 \ where C 
is replaced by C s and C s+1 , respectively. However, in this case the lower 
bounds 1/16 for A and 1/24 for C are to be replaced by 1/6 and 1/9, 
respectively. 

The closer A and C to optimal values 1/16 and 1/24 are, the smaller x 
is to be chosen. The inverse is also correct. Namely, in the next subsection 
we prove that estimates (2.9) and (2.10) become asymptotically exact in 
the limit x ~^ 0. I n this case factorials (3s)! and (3s + 3)! in the right- 
hand sides of (2.26) can be replaced by other expressions g(s) and h(s) that 
provide more precise estimates for R^ q \ Indeed, repeating the computations 
of subsections 2.3.3 and 2.3.4, one can see that in the limit x ~~ * function 
g(s) can be chosen close to (2s — 1)!!/4 S . This make an evidence for the 
central limit theorem to hold for the centered random variables 

NL° a = TrH a - E{ TrH a }. 

This observation explains also the fact that the odd "moments" of the vari- 
able L° decrease faster than the even ones as N — > oo. That is why the 
estimates for R( 2s ^ have the form different from those of _r( 2s+1 ) and are 
proved separately. 

For finite values of x, the use of some expression proportional to (3s)! is 
unavoidable. 

2.5 Proof of Theorem 2.2 

We present the proof of Theorem 2.2 for the case when k is fixed and N — > oo. 
Regarding relation (2.19) with q = 2, we obtain relation 

^ - \ K> • M) 2fc _ 2 + ^ • 2 -^>M 2t _, + i D g>_ 2 . ,,44) 

Proposition 2.1 implies that = 0(N~ 4 ) and that M^f ) -m k = 0(N- 2 ). 
Then we easily arrive at the conclusion that 

Agf-^ + O^). (,45) 

where r k are determined by relations ro = and 

1 x 1 2A;(2A; - 1) , 
r k = - ( r * m ) fc _ 1 + _ . L mk _ u k >\. (2.46) 
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Passing to the generating functions and using relations (2.11) and (5.11), 
we obtain equality 



1-r) 



k 



Returning to relation (2.8), we conclude that 
rW , 1 _(2) 



M. 



2k 



1 

N 4 



^ + Jf2 m k ^ I — ) ' 



Indeed, the difference between ML and rrik is of the order N and the 
next correction is of the order iV -4 . Regarding mj^ and using (2.45), we 
obtain equality 



m 



(2) 



1 

fc-i + 16 



{l-rf 



k > 1, 



(2.47) 



J fc 



and m^ 2 -* = 0. Solving (2.47) with the help of (2.11), we get expression 

^2 



m 



(2) 



1 

16 



(1 - r) 5 /2 

It is easy to see that (2.48) implies relation 

1 (2fc-3)(2ife-2)(2A;-l) 



(2.48) 



1 
16 



(1 - r) 5 /2 



16 



3! 



-m k -2 



and hence (2.6b). Theorem 2.2 is proved. 



2.6 More about asymptotic expansions 

The system (2.22)-(2.23) of recurrent relations is the main technical tool in 
the proof of the Proposition 2.1, where the estimates for B and R are given. 
However, the crucial question is to find the correct form of these estimates. 
The first terms of the asymptotic expansions described in previous subsec- 
tion give a solution of this problem. Indeed, repeating the proof of Theorem 
2.2, we see that formulas (2.46) and (2.48) indicate the form of the estimates 
to be proved. Then the proof of Proposition 2.1 is reduced to elementary 
computations, where the most important part is related with the correct 
choice of the factorial terms in inequalities (2.15). 

The next observation is that relation (2.23) resembles inequality (2.20) 

obtained from (2.19) by considering the absolute values of variables D^... t a q 
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and replacing in the right-hand side of (2.19) the sign "— " by the sign "+". 
So, relation (2.23) determine the estimating terms with certain error. 
However, it is not difficult to deduce from estimates (2.25) and (2.26) that 
if q = 2s, then this error is of the order smaller than the order of R( 2s \ 
This means that relations (2.23) determine correctly the first terms of the 
1/iV-expansions of all R( 2s \s > 1 and not only of R^ as mentioned by 
Theorem 2.2. The same is true for the 1/N expansions of D^*\ It is easy 
to show by using (2.23) and results of Proposition 2.1 that these corrections 
are given by formulas 

where are such that the corresponding generating function f( 2s )(r) = 
J2k>o r^T h verifies equation 

r^(r) = ^f<*>(r) + (2, - l)r^{r) (r f \r)) . (2.49) 

Using equalities (2.11) and resolving (2.49), we obtain expression 

(2s) (2,-1)!! 

k ( 4A T2)s 



(1-T 



,2s 



The left-hand side of relation (2.23) for R^ involves variables Rj 9 \ 

Rj Q 1 \ and R^ q+1 \ This can lead one to the idea to use the generating 
functions of two variables G(r, fi) to describe the family of numbers R. In 
this connection, the following comment on the structure of the variables 
could be useful. Introducing a generating function F(t) = X^>o rJ A?; we 
see that 

J2D$T 2k = E{[F°(T)}<>}, 

k>i 

where F°(t) = F(t) — EF(r). Then the mentioned above function can have 
the form 

G D {r^)= £ Z^ 2 ^ = E{e^)}-l. 

In particular, regarding such a generating function of r^ s \ one arrives at 
the expression 



r> t \ ( 2s ) 2k M 2s j A 1 

G r (r,n)= ^ r\W =exp - 



fc>l,s>l 



(2s)! P \4iV2(l-T 



2 
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This expression show that the central limit theorem can be proved for the 
random variable NF°(t) in the asymptotic regime k 3 /N 2 <C 1 mentioned in 
Theorem 2.2. This asymptotic regime can be compared with the mesoscopic 
regime for the resolvent of and the central limit theorem valid there [2] . 



3 Orthogonal and anti-symmetric ensembles 

In this section we return to Hermitian random matrix ensembles H^) with 
rj = 1 and r/ = — 1 introduced in section 2. Let us consider the moments 
of H^> . Using the method developed in section 2, we prove the following 
statements. 



Theorem 3.1 (GOE). 

Given A > 1/2, there exists \ such that 



MW<m k + A± CU) 



for all k,N such that k < k^ and (2.13) hold. If k is fixed and N — > oo, 
then 

=m k+ l -^ k +J)^ + 0(A r-i) (3.2) 



2k ■ 2N 

and 



D% N) = £ E{L° a L° b } = ^ + 0(N-*). (3.3) 



a+b=2k 



The proof of Theorem 3.1 is obtained by using the method described in 
section 2. Briefly saying, we derive recurrent inequalities for and D^ k \ 

then introduce related auxiliary numbers B and R determined by a system 
of recurrent relations. Using the triangular scheme of recurrent estimates 
to prove the estimates we need. Corresponding computations are somehow 
different from those of section 2. We describe this difference below (see 
subsection 3.1). 

Let us turn to the ensemble H^ _1 \ Regarding the recurrent relations for 
the moments of these matrices, we will see that for M 2 ( £" 1} are bounded 
by M^T 1 ^ . Slightly modifying the computations performed in the proof of 
Theorem 3.1, one can prove the following result. 

Theorem 3.2 (Gaussian anti-symmetric Hermitian matrices). 
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Given A > 1/2, there exists x > such that the moments of Gaussian 
skew-symmetric Hermitian ensemble ^ admit the estimate 



M^><m k + A- (3.4) 
for all values of k,N such that (2.13) holds. Also 

\ D *\ = °(w) and * i= °(f)' (3 - 5) 

Given k fixed, the following asymptotic expansions are true for the moments 
ofH(-V 

i\A N ) , 4,0 - (fe + IK , , AT -u fo R \ 
M K 2k ' = m k + — +o(N ), (3.6) 

and for the covariance terms 

D% N) = £ E{^ i L° 2 } = ^±l + 0(iV- 3 ). (3.7) 

ii+C2=2fc 



3.1 Proof of Theorem 3.1 

Using the integration by parts formula (5.7) with n = 1 and repeating com- 
putations of the previous section, we obtain recurrent relation for M 2k = 
EL 2k ; 

1 2k — 1 1 

M 2k = - £ M 2k _ 2 _ 2j M 2j + — — .U 2 , 2 + t £ E K^U • 

j=0 ai+a 2 =2fc-2 

(3.8) 

Regarding the variables 

2fc 2fc 

£ ^S!...^= E v{L° ai L° a2 ... L ° aq } 

ai,...,a q ai,...,a q 

and using formulas (5.6) and (5.8) with rj = 1, we obtain relation 

ai—2 1 ai— 2 



r>(<?) — i V A/f n^) j- _ /t/f 

-^ai,...,a„ " 4 L lvl 3 ±J a 1 -2-j,a 2 ,...,a q t 4 L -Z-J^jw ,...,a, 

j=0 j=0 

^ ai —2 ai —2 

+ 7 E £, jS-2-j,o 2l ..,o,-7 E ^i 2 ii-2-i^iLl ) a 9 + 7^( ai ~ 1 ) E { L a 1 -2 L a 2 
J"=0 j=0 
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1(1) 

^2 a i^ai+a i -2-Da2,...,a i _ 1 ,a l+ i,...,a 9 + r ) ^2 X] ai ^d2,-,«i-i ,a;+ai -2,o i+ i ,...,a g • 

(3.9) 



2_/V 2 r - "' "it"! - z ^a2,...,aj_i,a» + i,...,aq 1 27V 2 



Introducing variables 

2fc 



£2= E 



oi,...,a q 

we derive from (3.9) inequality 

^ < \ e ^t 2 - 2 ,M 2j + ^ri + J e ^ri-2,^ + 4^- 2 

j'=o j=o 
g-l^.( ? -2) (2j + 2)(2j + l) (g-l)fc 2 

j=0 

We have used here the same transformations as it was used when passing 
from equality (2.19) to inequality (2.20). 

Now we proceed as in Section 2 and introduce the auxiliary numbers B 
and R that verify relations 

B k = ±{B*B) k _ 1 + ±B k - 1 + ±RW v k>l, (3.11) 



k-i-j 



and 

4«) = i( B , fl «) t _ i + ti( B »,^» 

The initial conditions are: B = 1, i^ 2) = l/(2iV 2 ). The triangular scheme 
of recurrent estimates implies inequalities 

M^<Bf\ and \D$\ < < . (3.13) 
The main technical result for GOE is given by the following proposition. 

Proposition 3.1 

Let us consider B and R for the case of GOE (rj = I). Given A > 1/2 and 
1/4 < C < 2 • 6! ; there exists x such that the following estimates 

B ( k N) <m k + ^, k>2, 
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or equivalently 



B { k Q) < [f(r)] k + ^ 



A 

N 



l-T 



, k>2, 



and 



and 



(2s) C(3s)! 

U k S N 2s 



(2s+i) C(3s + 3)\ 

tt k - N 2s+2 



(l-T) 



2 s 



(1 _ r )2s+5/2 



(3.14) 

(3.15a) 
(3.156) 



hold for all values of k,q and N such that k < and (2.13) and (2.16) 
hold. 



The proof of this proposition resembles very much that of the Proposition 
2.1. However, there is a difference in the formulas that leads to somewhat 
different condition on A. To show this, let us consider the estimate for B k . 
Substituting (3.14) and (3.15) into the right-hand side of (3.11) and using 
(2.11), we arrive at the following inequality (cf. (2.28)) 



A 
N 



k Ak 



l-T 



+ 



A 2 + 6C 
4N 2 



(l-T) 



that is sufficient for the estimate (3.14) to be true. Taking into account that 

= km k -i, (3.16) 



we obtain inequality 



. 1 2A + A 2 + 6C 
A> — | . 

~ 2 4iVm fc _i 



It is easy to show that m^^\\fk > (2k) 1 . Then the last inequality is 
reduced to the condition 



A> - + (A + A 2 + 3C)^l. 



(3.17) 



The estimates for also include the values y/x and %. We do not present 
these computations. 
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Let us prove the second part of Theorem 3.1. Regarding relation (3.8) 
and taking into account estimate (3.15a) with q = 2, we conclude that 



M. 



(JV) 



m k + —m£ ) +o{N 1 ), as N -► oo. 



2k 



It is easy to see that the numbers mj^ are determined by relations 



m 



(i) 



1 / m \ 2fe - 1 

2l m * m J fc -i + ^— m *"i 



(3.19) 



and mg 1 ' = 0. Passing to the generating functions, we deduce from (3.19) 
equality 



m 



(i) _ 



r 


1 


1 - Vl - r 


1 (k + l)mfc 


.1-T. 


fc~ 2 


VI -T 


T 2 2 



Relation (3.2) is proved. 

Let us consider the covariance term D^> . It follows from the results of 
Proposition 3.1 that 

Then we deduce from (3.12) with q = 2 that r k is determined by the following 
recurrent relations 

1. , 12ife(2A;-l) 
r fc = - (r * m) fc _! + m fc _i. 

Solving this equation, we get 



-2(1 -r)2 

This completes the proof of Theorem 3.1. 



3.2 Proof of Theorem 3.2 

In present section we consider the ensemble with rj ■ 
the elements of H (2.3) are given by imaginary numbers 



T. In this case 



(H) X y - 



--W x 



x < y 



Jxy ~ s/N V ' 
and the skew-symmetric condition holds: 

{H) X y = —{H) yx = {—H) yx 
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Regarding the last term of the formula (5.7) and using previous identity, we 
can write that 



2k-l N -. 2fe-i ( 1 s 

4^ E E E{(W-% x (H 2k -^) yx } = ~\ £ (-ir^j^ Tr^-n 

i=l x,y=l j=l 



Then we derive from (5.7) equality 

E^ 2 fc = 4 E E{LjL 2 k-2-j} ~ —~EL 2 k-2 

3=0 

that gives recurrent relation for the moments of matrices ^ 

M 2k = - ]T M 2j M 2jk _ 2 _ 2j - ^M 2fc _ 2 + Z D^_ 2 , (3.20) 

j=0 

where the term is determined as usual. Regarding the general case of 
D(i\q> 2, and using (5.8), we obtain the following relation 

m-2 / i \ai +1 



D ai,...,a q = 4 E E {- L ai-2-j-t'j[^a 2 " " " L aJ°/H ^ E { L ai-2 I 'o 2 " " " L a 9 

+ 1+ 4 ^T +1 E E { L a 2 ---^_ 1 ai W ai _ 2 ^ +1 ...L° 9 }. (3.21) 

i=2 

Comparing this equality with (3.9) and then with (3.10), we see thatM 2 ( J?" 1} 

and D^ k ,ri ~ ^ are bounded by the elements B and of recurrent rela- 
tions determined by equalities (3.11) and (3.12), respectively. Indeed, tak- 
ing into account positivity of M 2 j and the fact that I /AN < k/4N and 
1 + (— l) a » +1 < 2, we obtain inequalities 



M 2k < l -Y, M 2j M 2k _ 2 _ 2j + ^M 2fc _ 2 + \D 

3=0 



(2;u=-i) 

2fc-2 > 



where 



= (q tV =-l) , j=.(q,r)=l) 
U 2k - U 2k 



Then we conclude that 



Mt _1) < Mt 1] < B k and < J#>. (3.22) 
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This completes the proof of the first part of Theorem 3.2. 



(2) 

Now let us turn to the asymptotic expansion of M 2 k and DV for fixed k 
and iV — > oo. Regarding (3.21) with q = 2, we obtain the following relation 



a\+a,2=2k 



AN 2 

Now, introducing variable r k 



, 1 + (-!) a2+1 V nVT 

ai+a2=2fc 



and taking into account estimates (3.4) and (3.5), we conclude after simple 
computations that r is determined by recurrent relations 

r k = -( r *m) k _ 1 + —m k - 1 . 



Using relation (5.10), we can write that 

k 2 2k(2k - 1) k 1 

— m k -\ = — m k -i + -m k -i = - 



2 

Then 



2-4 



(1 - r) 3 / 2 



1 

+ 4 



{l-rf 



1 

+ 4 



1-T 



(l-r)V2_ 

(3.23) 



Now let us consider l/A^-expansion for M 2k 

M 2 k = m k + j^m^K 



It is easy to see that equality (3.20) together with estimates (3.5) implies 
the following recurrent relation for m k 



(i) 



m 



fc 2 V >k-\ 



m k -i 



Then 



m 



(i) 



rf(r) 

4 v / r^ 



4,0 ~(k + l)m k 



(3.24) 



31 



These computations prove the second part of Theorem 3.2. 



Theorems 3.1 and 3.2 show that there exists essential difference between 
GOE and Gaussian anti-symmetric ensemble. Let us illustrate this by the 
direct computation of for these two ensembles. 

In the case of GOE, we have 

1 V Fff 2 = lVFl/ 2 + IVFy 2 _ N(N-1) 111 
^n xy n2 ^ n v xy t n2 2_ ^ v xx ^ N 2 ^ 2N 4 4N' 

x,y=l x<y x=l 

This relation reproduce (3.2) with k = 1. 

The first nontrivial moment of anti-symmetric matrices reads as 

iVEtf = lVE^ = N(N-1) = 1__J_ 
N Z - 1 xy N 2 2-? x v 4iV 2 4 47V 

x,y=l x<y 

that agrees with (3.24). 

Finally, let us point out the difference between GOE and anti-symmetric 
ensemble with respect to the first term of the expansion of D( 2 ) given by 
(3.3) and (3.23), respectively. This indicates that Gaussian Hermitian anti- 
symmetric ensemble represent a different universality class of the spectral 
properties of random matrices than that of GOE (see for example, the mono- 
graph [16]). 



4 Gaussian Band random matrices 



Now let us consider the ensemble of Hermitian random matrices given by 
the formula 

#(iv,&) 

where {h xy , x < y} are determined by (2.1) and (2.2) with r] = 0. The 
elements of non-random matrix U = U^ N ' b ^ are determined by relation 




x,y = l,...,N, 



where u(t), t £ R is a positive even piece-wise continuous function such as 

/oo 
u(t) dt = u\. 
-oo 



Til 



x,y = l,...,N, 



(4.1) 
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Without loss of generality, we can consider uq = 1. We assume also that 
u(t),t > is monotone. If u(t) is given by the indicator function of the 
interval (—1/2, 1/2), then matrices (4.1) are of the band form. We keep the 
term of band random matrices when regarding the ensemble (4.1) in the 
general case. 

It is known (see for instance [5, 17]) that the moments of H^ N ' b ^ converge 
in the limit of 1 <C b <C N to the moments of the semicircle law; 

m k (ui), (4.3) 

where the numbers {mk(ui), k > 0} are given by recurrent relations 

fc-i 

mk(ui) = — m k -i-j{ui) mj(«i), m (iti) = 1. 

3=0 

The generating function /i(r) = T>t m k (ui ) is related with /(r) (1.7) by 
equality /i(r) = /(rui) and therefore 

m k (ui) = u\m k . 



M, 



(N,b) 
2k+l 



0, M 



2fc 



4.1 Main results 

In this section we present non-asymptotic estimate for the moments of 
M 2 ( f' &) . This improves proposition (4.3). Let us denote 



-(b) 1 . 1 V 2 ? 

(=— oo 



Clearly u\ > u\ and — ► u\ as 6 — > oo. 
Theorem 4.1. 

Given a > 1/12, i/iere existe 6* > stic/i that the estimate 

M^ b) <(l + cm^^-J mfcCth), (4.4) 

(fc + l) 3 

where u = max{ui, 1/8}, ZioWs /or a// values of k, b such that ± — ^ < 6 
and b< N. 
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The proof of this theorem is obtained by using the method described 
in section 2. We consider mathematical expectations of variables L 2k (x) = 
(H 2k ) xx and derive recurrent relations for them and related covariance vari- 
ables. Certainly, these relations are of more complicated structure than 
those derived for GUE in section 2. However, regarding the estimates for 
M 2k = EZ/2fc(x) by auxiliary numbers B k , one can observe that equalities for 
B k and related numbers is almost the same as the system (2.22)-(2.23) 
derived for GUE. This allows us to say that the system (2.22)-(2.23) plays 
an important role in random matrix theory and is of somewhat canonical 
character. The estimates for the moments M^' b ^ follow immediately. 

4.2 Moment relations and estimates 

In what follows, we omit superscripts (N, b) when no confusion can arise. It 
follows from integration by parts formula (5.8) that 

E{H xy {H l ) yx } = --U xy £E{(fF%, {H l -i) xx } . (4.5) 

j=0 

Then, regarding L k {x) = (H k ) xx , we obtain equality 

j 2fc-2 

EL 2fc (x) = E{L 2k - 2 -j(x) Lj[x}} , 

3=0 

where we denoted 

, N / _ \ 
y=l 

Introducing variables M k {x) = ELfc(x) and M k [x] = EL/Jx], we obtain 
equality 

k— 1 

M 2k {x) = - £ M 2k _ 2 _ 2j {x) M 2j [x] + -D$>_ 2 (x, [x]), (4.6) 

3=0 

where we denoted 

D$>_ 2 {x, [x]) = ]T E K» KM ■ (4-7) 

ai+a2=2k— 2 

In (4.6) we have used obvious equality M 2k+ i(x) = 0. 
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To get the estimates to the terms standing on the right-hand sides of 
(4.6) and (4.7), we need to consider more general expressions than M and 
D introduced above. Let us consider the following variables 

i<-*>(x) = E {(l/^flw • • • % r H^) x J , (4.8) 

where we denoted ir r = (pi,P2, ■ ■ ■ ,Pr+i) with Y%=i Pi = the vector 
Vr = (yi, ■ ■ ■ ,Ur) and ^> y denotes the diagonal matrix 

= S st U (t^j , s,t = l,...,N. 

One can associate the right-hand side of (4.8) with 2k white balls separated 
into r + 1 groups by r black balls. 
The second variable we need is 

D^ a tit\x q ) = E{ L° ai ( Xl )L° a2 [x 2 ] ...L° aq [x q ] }, (4.9) 

S v ' 

where a q = (a±, . . . ,a q ) and x q = (x±, . . . , x q ). We also denote \a q \ = 
J2f=i a i- So, we have the set of \a q \ white balls separated into q boxes by 
q — 1 walls. 

The brace under the last product means that the set {ai|a2| • • • \a q } of 
walls and white balls is separated into r + 1 groups by r black balls. The 
places where the black balls are inserted depend on the vector a q . 

Let use derive recurrent relations for (4.8) and (4.9). These relations 
resemble very much those obtained in section 2. First, we write identity 

N 

M^\x) = £ E {H xs (H^^ yi H^ ■ ■ ■ % r H p ^)J , 

s=l 

and apply the integration by parts formula (4.5). We obtain equality 
M£'*\x)= £ B{L ai (x)L a2 [x]}. 

ai +a 2 =2k-2 s ~y ; ' 

In this relation the partition ir' is different from the original it from the left- 
hand side. It is not difficult to see that it' depends on particular values of a\ 
and a2, i.e. on the vector (ai, 02). Returning to the denotation M = E{L}, 
we obtain the first main relation 

M^ r \x)= £ M ai (x)M a2 [x}+ J2 D$$»<»)\x,[x]). 

ai+a 2 =2fc-2 V ~^ ! ' ai+a 2 =2fc-2 

(4.10) 
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Let us consider 



N 



sxi 



=1 ^- 



L° a2 [x 2 ]---L°[x q ] } 



and apply (4.5) to the last mathematical expectation. We get 

i 0,1— 2 

D^ltiV^) = 7 E E { L ai _ 2 _ a ,(xi)L a ,[xi] \K 2 [x 2 }...L aq [x q ]\° } 



a'=0 



<(yr,a' q+1 ) 



+7I2E E E { (^*x i fl tt *- 1 -^x 1 ^- 1 ) seixi ^N- Xi---K q [ Xq ] }. 



i=2 j=0 



< +2 (?/r + 2><'+lW) 



In these expressions, 7r' and it" designate partitions different from 7r; they 
depend on the vectors a! q+1 = (a\ — 2 — a', a 1 , a 2 , ... , and 

"g+iOO = C?> a i - 1 - j,ai - l,a 2 , . . . ,ai_i,aj_i, . . . ,a 9 ), 

respectively; also y^.+2 = (a^i , a^i , 2/i , J/2 5 • • • ,Vr)- The notation Xj in the last 
product of (4.11) means that the factor L ai is absent there. Repeating the 
computations of section 2, we arrive at the second main relation 



f)(q,Tv r (a q ),y r )(- \ _ j, 



6 



(4.12) 



where 



ai-2 



T l = 7 E ^ai-2-a'(^l) ^t 2 ,..,a,([^]'N,...,N) 



a'=0 



ai-2 



<(j/r,a' 9+1 ) 



T 2 = ^ E • • • , [Xj); 



o'=0 



1 v^ 2 n (9+ 1 ><(j/r-,a; +1 )) . r i r i r n 

T 3 = 7 E £ 'a 1 -2-a',a',a 2 ,...,a 9 ( a; l ) M> M, • • • , [x,])j 



<(j/r,a' g+1 ) 



a'=0 
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a'=0 s ~v ' 

T 5 = 712 E E M ai+ai _ 2 (x!) 4^-1^+1,..,^ (M, • • • , [Xi-l], [x i+ i], . . . , [s,]); 
40 i=2 i=0 v ' 

<' +2 (5r + 2' Q! ','+2( i )) 

and finally 

T 6 = 712 E E ^i!+ar-Tatt^iS+i,..,a 9 ( X l' H, • • • , ^i-l], • • • , 

i=2 j=0 

Now let us introduce auxiliary numbers {B^ N ' b \ k > 0} and 

^* ;JV,6) = ^Kfi, for g>0anda,>0, 

determined for all integer k, q and a» by the following recurrent relations (in 
B and .R, we omit superscripts N and b). Regarding {B}, we set Bq = 1 
and determine .B^ by relation 

^ = tEVi-A + 7 E ^Sa 2 , *>1- (4-13) 

j=0 ai+a 2 =2fc-2 

Regarding {R}, we set i?^ ^ = 1 and = 0. We also assume that 
Ra] = when either q > \a q \ or one of the variables Oi is equal to zero. The 
recurrent relation for R is 

* ai— 2— j 

-. ai—2—j ai—2—j 

+i V +- V i? (2) 

T 4 JL 7,ai-2-j,a2,...,a 9 ~ 4 Jl ai-2-jj"- I, '02,...,o, 

^4i,2 Zw J £ 'ai+a/-2- n 'a2,...,ai_i,Oi + i,...,o, ^ 4 l2 "*- n, 02,...,ai_i,Oi + i,...,a 9 ■ 
i=2 i=2 

(4.14) 

Existence and uniqueness of the numbers B and R follow from the triangular 
scheme described above in section 2. 
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Using the triangular scheme of section 2, it is easy to deduce from rela- 
tions (4.10) and (4.11) that 



supM 2 ( £ r,Sr) (x) < B k 



and 



sup \D$ZfcM (xi, [x 2 ], . . . , [x q ]) | < R$ a2 _ 

Xq ,y r 



(Iq ' 



(4.15) 



(4.16) 



Let us note that when regarding (4.15) with k = 0, we have used the property 
of u (4.2) 



u 



i=l 



x-Vi 



<u r <l. 



Now, let us introduce two more auxiliary sets of numbers B k and R(q) k - 
We determine them by relations 



^ = 7"E B k-i-jBj + jRk-n B = 1, 

3=0 



fc-1 



1 



(4.17) 



and 



4" = | E + ^ £ (2i + 2>(2i + 1) • 

j=0 j=0 
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It is clear that 



B k <B k and ]T tfg...^ < R 

ai+ . . .a q =2k 

The main technical result of this section is as follows. 



(?) 



(4.18) 



(4.19) 



Proposition 4.1. 

Lei « = max{tti, 1/8}. Given A > 1/16, there exists 9 > swc/i £/iat t/ie 



h(r) + 



Aii 



b 2 (l-rni) 5 /2 



(4.20) 



ZioWs /or a// values of k < ko, where ko verifies condition k^ < 9b 2 . Also 
there exists C 

1 3^4 
— <C<max{ — ,4!} 
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such that inequalities 



fc 1,2s 



(1 - r«i) 2s 



and 



(2s+1) ^+1(35 + 3)! 

- ° 



(1 - TUl) 2 ^ 1 



(4.21a) 
(4.216) 



hold for all values of k and s such that 

2k + q< 2k 
with q = 2s and q = 2s + 1, respectively. 

The proof of this proposition can be obtained by repeating the proof of 
Proposition 2.1 with obvious changes. The only difference is related with the 
presence of the factors u\ in (4.17) and (4.18). This implies corresponding 
changes in the generating functions used in estimates (4.20) and (4.21). Also, 
the conditions for A (2.29) and C (2.38), (2.43) are replaced by conditions 

. 3C A 2 Ul n 
C > 6sA + C(1 ~ 5sa) +26(1 + 10nC(l + C) + 2u\AC) 



and 
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179m > 20 + 3nC + IWuuiA. 



The last inequality forces us to use u instead of u\ in the proof. Otherwise, 
we should assume that u\ > 1/8. We believe this condition is technical and 
can be avoided. 



4.3 Spectral norm of band random matrices 



Using this result, we can estimate the lower bound for b to have the spectral 
norm of \\H^\ 



Xilmx bounded. 



Theorem 4.2 If 1 < (logiV) 3 / 2 < b, then A 
1. 



(N,b) 
max 



ui with probability 



Proof. Using the standard inequality 
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we deduce from (4.4) estimate 

A 4k + l? \ k 
P >V^i(l + e)}< A V ^ (1+g 6 )2fc J 4 (4.22) 

that holds for all k + 1 < fl 1 / 3 ^ 2 / 3 , where is as in Theorem 4.1. In (4.22), 
we have used inequalities m^{u\) < 7«2fc and rri2k < 1- 

Assuming that 6 = </>Ar(log A") 3 / 2 , where ^Tv^ooasA^oo, and taking 
k + 1 = t6» 1/3 6 2/3 , < t < 1, we obtain the estimate 

^{A r ( ^ x 6) > v^TU + e)} < Aexp{-2^ 1 / 3 & 2 / 3 log(l + e) + 2ant 3 }-^^. 

(4.23) 

Using relation ui = ui(l + 1/6), we easily deduce from (4.23) that 

^{A^ak 6) > v^I(l + £)} < N 1 -^ 1 ^ 3 

with some positive C. Then corresponding series of probability converges 
and the Borel-Cantelli lemma implies convergence of X^ax to y/ui. Theorem 
4.2 is proved. 

Let us complete this subsection with the following remark. If one opti- 
mizes the right-hand side of (4.23), one can see that the choice of t = to = 
6 1 / 3 v /log(l + e)(au)~ l l 2 6~ l l ? ' gives the best possible estimate in the form 

7Vexp{-&^-L=(log(l + e)f' 2 }. 
v2au 

Once this estimate shown, convergence — > \Ju{ would be true under 
condition that b = O(logA). However, one cannot use the optimal value 
of to mentioned above because this choice makes k to be k = 0(b). This 
asymptotic regime is out of reach for the method of this paper. 

5 Auxiliary relations 

5.1 Integration by parts for complex random variables 

Let us consider matrices H with elements H xy = v xy + iw xy , where the 
family {v xy ,w xy , 1 < x < y < A} is given by jointly independent Gaussian 
random variables with zero mean value. We denote 
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Let us assume that x < y. Then integration by parts formula says that 

{ dv xy } + l( " xyE { dw. 



It is easy to see that 
d(H%t _ >A A j-i dH s f t ' i-j TTj-i rri-j ,TTj-i Tji-i 

—fa, ~2s 2s H ss> • - R vt ~2s [ H sx H yt + H ty H xt 

ov *y j= i s >,t>=i OVx v j=i 



Similarly 

d(H% 



(5.2) 



= ^'Hi,.H ! ;ll J Hi/H 1 ,, 1 ■ (5-3) 



OW X y j = 1 

Substituting (5.2) and (5.3) into (5.1), we get equality 

i 

BH xy (H% = (^ - Cxy) E V{Hlx l Kt 3 } 

i=i 

i 

• • wE 1 '^ 1 '''/}- ,• y . (5.4) 

It is not hard to check that the same relation is true when x > y. Also 

l 

EH xx (H l ) st = £ xx ]T V\Ili,. Il',, >}. (5.5) 
i=i 

5.1.1 Gaussian Ensembles {H^)} 
Regarding formulas (2.1)-(2.3), we see that 

Vxy _ W^xy 

Vxv ~7n' Wxv ~7n 



and 



Regarding the sum of (5.5) with doubled (5.4), we obtain relation valid for 
all values of x and y 

EH xy (H l ) st = ^ E '//:,'} + ^ E HUH/"',! 1 }- (5-6) 
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Let us mention two useful formulas that follow from (5.6); these are 

1 1 

E Tr(H l+1 ) = — E { TrF^ 1 TrH l ~ j } 

3=1 

I N 

+i^E E E{(Hi-X x (H l -i) yx } (5.7a) 

j=l x,y=l 

and 

TrF< = E ji^][>') ss j = ^E^ 1 + ^H l ~\ (5.7b) 
5.1.2 Band Random Matrices 

Using (5.4) and (5.5) in the case of matrices (4.1), we see that 

t _ 1 + $xy TT > _ 1 ~ &xy TT 

Zxy ~ o u xy, <^xy — o u xy 



Then (5.4) and (5.5) imply equality 

EH xy (H% = KHIi/ll'-r 1 }- (5.8) 

i=i 

Regarding this relation, one can easily obtain analogs of formulas (5.7a) and 
(5.7b). 

5.2 Derivation of Equality (2.18) 

We consider the case of Hermitian matrices rj = only. Regarding (2.17), 
we can write that 

E{L° ai ..-L° aq } = E{L ai Q°}, 

where Q = L° 2 • • • L° . Using integration by parts formula, we obtain as in 
(5.1) that 

E{H xy (H^) yx Q°} = k»E j + i^E | } • 

Obviously, 

y x V = V E I LfP H ai ~ 1 - j + W H^' 1 -^ Q°\ + H 111 - 1 —^- 
dv xv f-f I L + ra ^ J V J + yx dv xv ' 
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It is clear that 



and 



dV X y dV X y 



OV X y 1\ 



iV 



fja.i—1 I frai-1 

xy ' yx 



Also we have 



dm^Q 



ai — 1 



UW X y J UUJ X y 



It is clear that 



and 



8W X y 8W X y . =2 



W "' / ' ; '' 



-= *jEE >^ ; /0 ' - "i;'"T J 



dw xy N t=lj=1 



N 



jjai — l _ jjai — l 



Gathering these terms, we finally obtain that 

E{H xy (H^-% x Q°} = -L a f^v{H^H£-i-iQ°} 



3=1 



1 q 

4j^E aiE { ff P lL °a 2 - ■ ■ K^Hxy 1 L °a l+1 ■ ■ ■ L a q } ■ 
i=2 



Now (2.18) easily follows. 

5.3 Catalan numbers and related identities 

In the proofs, we have used the following identity for any integer r > 1, 

(2k+2r\ 



(1 - r) r +V2 



^2k+2r^ 
r (k+r\ mk ' 



(k+r\ 
\k+V 



(5.9) 
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or in equivalent form, 



1 



1 (2k + 2r)\ 



r\ 



_(l_ r )r+l/2 

Two particular cases are important: 
(2A: + 2)(2fc + 1) 



2 2k k\ (2r)! (k + r)V 



and 



(2/c + l)(2/c + 2)(2fc + 3) 



3! 



We also use equality 



1 



(1-r) 



l+i 



(1-r) 3 / 2 . 
1 

(l-r)5/2_ 
(A; +!)■■■(*; + Z) (fc + Z)! 



Z! 



fe! /! 



(5.10) 

(5.11) 
(5.12) 

(5.13) 
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